Renormalization of hole-hole interaction at decreasing Drude conductivity 
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The diffusion contribution of the hole-hole interaction to the conductivity is analyzed in gated 
GaAs/Ina;Gai_i,As/GaAs heterostructures. We show that the change of the interaction correction 
to the conductivity with the decreasing Drude conductivity results both from the compensation of 
the singlet and triplet channels and from the arising prefactor a, < 1 in the conventional expression 
for the interaction correction. 



The quantum corrections to the conductivity in disor- 
dered metals and doped semiconductors are intensively 
studied since 1980.— Two mechanisms lead to these cor- 
rections: (i) the interference of the electron waves prop- 
agating in opposite directions along closed paths (WL 
correction); (ii) electron-electron (e-e) or hole-hole (h-h) 
interaction. 

The role of the e-e (h-h) interaction has been a sub- 
ject of theoreticalii^ii^i^i^i^i^ii^i^ and experimentaUSiiiii^ 
studies for more than two decades. The new interest 
in the matter is associated with discussion of the nature 
of metallic-like temperature dependence of the conductiv- 
ity observed at low temperature in some two dimensional 
(2D) systems, e.g., in n-Si MOSFET and in dilute 2D hole 
gas in Al3;Gai_a;As/GaAs and Gei_j:Sia;/Ge structures 
(see Refs. [l3lT3fl5l and references therein). As a rule 
such behavior is observed in low-density high-mobility 
structures with the relatively large value of the gas pa- 
rameter rg ~ \/2/ {aBkp) characterizing the interaction 
strength and by too high value of Tr (hereafter we set 
h = 1, ks — 1), where as, kp, and t are the Bohr 
radius, the Fermi quasimomentum, and the transport re- 
laxation time, respectively. The role of the interaction at 
Ts > 3 — 5 (i.e., at strong interaction), and/or Tr > 1 
(i.e., at intermediate and ballistic regimes) was theoreti- 
cally studied in Refs.'5,'3'5','li|l3Iii[i9p9, the experimen- 
tal situation was reviewed in Refs. (I3|)l4j[l5|. 

It should be noted that the metallic-like behavior is 
observed when the conductivity is not too high, there- 
fore the corrections can lead to essential change of the 
conductivity with the temperature. The changing of the 
interaction correction at decreasing temperature and/or 
conductivity was theoretically studied in framework of 
the theory of the renormalization group (RG) in the 
papersi^iii^ii^iii It has been shown that the correction 
renormalization depends on both the Drude conductiv- 
ity and the Fermi liquid amplitude 72 that controls the 
e-e interaction in the triplet channel. The contributions 
from singlet and triplet channels are opposite in sign fa- 
voring localization and antilocalization, respectively. In 
conventional conductors with high values of the Drude 
conductivity, cto = nkplGo 3> Go [where I is the mean 
free path and Go — e^/(27r^?i) ], the initial value of the 
amplitude 72 is small, and the net effect is in favor of 



localization. At cto ^ (10 — 15)Go or in dilute systems, 
however, this amplitude may be enhanced due to e-e cor- 
relations and thus results in metallic sign of da/dTJ^^ 

Significantly less is known about the role of the in- 
teraction correction in disordered 2D systems when the 
kpl value tends to unity, i.e., at crossover from weak 
to strong localization. Experimentally, this effect was 
studied in the simplest single-valley electron 2D system 
GaAs/Ini_j;Gaa;As/GaAs with small g factor.^^ It was 
shown that the net value of the interaction correction de- 
creases rapidly with the (Tq decrease at ctq ^ (12 — 15) Gq 
{kpl ^4 — 5). Such a behavior can result from the com- 
pensation of the contributions of the singlet and triplet 
channels as well as from suppression of both contribu- 
tions with decreasing ao- It is impossible to separate 
these two effects in the systems with small value of the 
g factor. The situation changes drastically when dealing 
with a system with large enough g factor. In this case the 
magnetic field can be used as a tool allowing to control 
the ratio between the two different contributions because 
it strongly suppresses the triplet channel and leaves the 
singlet channel unchanged. As shown below the hole 2D 
gas in strained GaAs/Ina;Gai_j;As/GaAs structures is a 
suitable object to study the renormalization of the in- 
teraction quantum correction with the conductivity de- 
crease. In a previous paper, Ref. [23 , we have studied 
these structures at high Drude conductivity, ctq > 30 Gq. 

In this paper we report the results of experimental 
study of the evolution of the interaction correction to 
the conductivity in a p-type 2D system with decreasing 
Drude conductivity within the range from ~ 30 Go to 
~ 3 Go when the ballistic contribution of the h-h inter- 
action is small. Firstly, we will outline the procedures 
used for extracting the diffusion part of the interaction 
correction and the value of the Fermi liquid parameter 
Fq — —72/(1 + 72) from the dependences of p^x and 
Pxy on the temperature and magnetic field. Then, we 
will discuss the change of Fq with decreasing Drude con- 
ductivity. Finally, we will show that the reduction of the 
interaction correction with the decreasing Drude conduc- 
tivity results from both the compensation of the singlet 
and triplet channels and from the arising of a prefactor 
Qfi < 1 in the conventional expression for the interaction 
corrections^ 
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I. EXPERIMENT 



We have measured the temperature and magnetic 
field dependences of p^x Sind p^y in the heterostructures 
GaAs/Ina;Gai_j;As/GaAs grown by metal-organic vapor 
phase epitaxy on semiinsulating GaAs substrate. The 
lattice mismatch between Ina;Gai_a;As and GaAs results 
in biaxial compression of the quantum well. The struc- 
tures consist of a 250 nm-thick undoped GaAs buffer 
layer, carbon (5-layer, a 7 nm spacer of undoped GaAs, 
a 10 nm Ino.2Gao.8As well, a 7 nm spacer of undoped 
GaAs, a carbon 5-layer and 200 nm cap layer of undoped 
GaAs. The samples were mesa etched into standard Hall 
bars and then an Al gate electrode was deposited by 
thermal evaporation onto the cap layer through a mask. 
Varying the gate voltage Vg we were able to change the 
hole density p and mobility p. within the following ranges: 
p = (2.5 . . . 8.0) X 10" cm-^ p. = (1000 . . . 5700) cmVVs. 
Two Hall bars prepared from each of the waffles 3856 and 
3857 with close parameters were measured. 

The magnetic field dependences of pxx and pxy at 
T — 1.4 K at different gate voltages for one of the sam- 
ples investigated are presented in Fig.[TJ It is clearly seen 
that despite the very large difference in conductivity val- 
ues at i? = 0, the magnetoresistance (MR) curves Pxx{B) 
are very similar: the sharp negative MR at low mag- 
netic field, which results from suppression of the interfer- 
ence contribution to the conductivity, is followed by the 
parabolic- like MR caused by the interaction corrections^ 
Since our goal is to study the interaction correction 
let us briefly explain the method allowing us to extract 
it from the experimental data. Under our experimental 
conditions the parameter Tr is small enough (Tr < 0.1) 
and therefore the main contribution comes from the diffu- 
sion part of the interaction correction. The unique prop- 
erty of the diffusion part is that it contributes to axx 
but not to (7xy This fact opens a possibility to extract 
this correction reliably even when the correction value 
is small. The most straightforward way is to find such 
contribution to axx which is absent in (Txy We extract 



these contributions by making use of the structure of the 
components of the conductivity tensor (Jxx and axy As 
shown in Ref. [U the weak localization correction and the 
ballistic part of the interaction corrections are reduced to 
renormalization of the transport relaxation time and can 
be accounted for through the temperature and magnetic 



field dependence of the mobility. Thus ax 
be written as 



and 



a xy can 



(Txx{B,T) 

axy{B,T) 



epp{B, T) 
1 + p^{B,T)B^ 

epp'^{B,T)B 
l + p^{B,T)B^' 



+ 6a':',{B,T), (1) 



(2) 



where Sa^^{B, T) is the diffusion part of the interaction 
correction. If the Zeeman splitting is very small as com- 
pared with the temperature, 5a^]^ is magnetic field inde- 




FIG. 1: The magnetic field dependences of pxx (a) and pxy 
(b) at T = 1.4 K for different gate voltages, which are charac- 
terized by the following values of p, ctq, and a{T = 1.4 K): 8 x 
10" 59.6 6*0, and 56.9 Go (curves 1); 4.5 x lO" cm-^, 

9.9 6*0, and 6.8 Go (curves 2); 3.9 x lO" cm-^ 8.1 Go, and 
4.37 Go (curves 3); 3 x lO" cm-^ 3.9 Go, and 0.36 Go (curves 
4); 2.6x10" cm"^ 3.5 Go, and 0.027 Go (curves 5). Structure 
3856. For clarity, the curves in the panel (a) are separated in 
vertical direction by the value of 0.2. 



pendent. It has the form^ii^iiii 

5cy^x'i{T) 



Go 



1 + 3 1 



In (1 + Fq" 
^0 



In Tr, (3) 



where the first term in square brackets is the exchange or 
the Fock contribution while the second one is the Hartree 
contribution (the triplet channel). For the following, we 
enter here the prefactor ai which was absent in Refs. 0, 
i0i 

Thus, knowing the hole density p we can find /i(T, B) 
from experimental axy vs B dependences [with the help 
of Eq. ([2])] and then calculate the first term in Eq. ([T|). 
The difference between experimental value of axx and 
this term should give the diffusion part of the h-h cor- 
rection to the conductivity. This method allows us to 
find Sa!^!^{B,T) for relatively low ao, when the interfer- 
ence contribution to MR is not negligible up to the high 
magnetic field. 

In what follows we demonstrate how this method works 
considering the results obtained for one of the samples, 
fabricated on the basis of structure 3857. 

Let us start with the case of high Drude conductiv- 
ity, ao 2^ 30 Go (for the details of determination of ao 
see Ref. 1 2 51). First for each temperature we have inverted 
the resistivity tensor whose components measured exper- 
imentally and found the conductivity tensor components 
axy and axx [solid curves in Figs.^Ka) and [2Kb)]. Then, 
using the obtained axy vs B dependences we have found 
Pl{B) [shown in inset in Fig. Ufa)] and calculated the ex- 
perimental value of the first term in Eq. ([1]). Finally, 
subtracting the latter term from the experimental value 
of axx we obtain Saxx [see Fig. djc)], which is identified 
with the diffusion part of the h-h correction da^^{B,T). 
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FIG. 2: (a) The experimental magnetic field dependences of 
a^y for two temperatures. The inset shows the magnetic field 
dependence of the mobility calculated from axy{B) for two 
temperatures with p = 5.4 x 10^^ cm~^. (b) The a^x vs B 
dependences for two temperatures. The solid curves are the 
data, the dotted curves are the first term of Eq. ([1} calculated 
as described in text, (c) The magnetic field dependences of 
the difference between experimental and calculated a^x val- 
ues for different temperatures, (d) The temperature depen- 
dences of /S.axx at different magnetic fields. The symbols are 
the experimental results; curves are calculated dependences 
with Fq = -0.4 (solid curves) and with Fq = -0.35 (dotted 
curves). Structure 3857, Vg = 2.4 V, ao ~ 30 Go. 
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FIG. 3: (a) The magnetic field dependences of a^x for two 
temperatures. The solid curves are the data, the dotted 
curves are the first term of Eq. calculated as described in 
text, (b) - (e) The temperature dependences of Acr^a; at dif- 
ferent magnetic fields. The symbols are the experimental re- 
sults; the curves are theoretical dependences calculated with: 
ai = 0.5, FS = -0.43 [panel (b)]; Ui = 1, F^ = -0.45 [panel 
(c)]; a. = 1, FS = -0.5 [panel (d)]; a, = 1, Fq" = -0.55 
[panel (e)]. Structure 3857, Vg = 2.8 V, p = 4.4 x lO" cm-^ 
o-o ==; 11 Go. 



these formulas, i o^^i^^ 



As seen from Fig.[5l^b) 5axx is a small difference between 
two large quantities. That is why an accuracy in deter- 
mination of 5a^l^{B ,T), i.e., the absolute value of the 
interaction correction, is sufhciently low. In particular, 
it is very sensitive to the value of hole density, which is 
experimentally known with some accuracy. However, the 
difference of the quantities 5axx taken at two tempera- 
tures for a given magnetic field (or taken at two magnetic 
fields for a given temperature) depends only slightly on 
the hole density and, therefore, is found with better ac- 
curacy. 

In Fig. [2Jd) we present the temperature dependences 
A(Ta;^(r, B) = 5axx{T, B) - 5axx{TQ, B), where Tq is the 
lowest temperature, obtained for different magnetic field. 
One can see that the higher is the magnetic field, the 
stronger is the change of Aaxx with the temperature. 
This dependence can be attributed to the Zeeman split- 
ting which leads to suppression of the triplet channel and, 
hence, to appearance of the magnetic field dependence of 
the interaction correction. Theoretically, the effect of 
Zeeman splitting has been considered in Refs. [8l l26ll27l . 
and 0. However, the expressions derived there are too 
complicated and, therefore, inconvenient for the practical 
use. Much simpler expression, which well approximates 
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^ + F-) 



J- Pi 



lnTr + 21nTrJl+(^^y 



(4) 



In Fig. m^d) we plot the curves calculated according to 
Eq. (gl) with = 1, g = 3^ and different values. 
One can see that the curves calculated with = —0.4 
almost coincide with the experimental data. 

Similar data treatment was carried out for the lower 
conductivity. In Fig. [3] we present the experimental and 
calculated magnetic field dependences of axx [Fig- E^a)] 
and Acr^jj-versus-T dependences for different magnetic 
fields [Fig. [3Jb)] for ctq = HGo {Vg = 2.8 V). As seen 
from Figs.[3l^c) ~[He), it is impossible to describe the data 
by Eq. ([4]) with the prefactor Ui = 1 for any F(f -values. 
This is not surprising because the theory predicts ai = 1 
only for large ctq value. However one can fit the data 
perfectly with a, = 0.5 and F^ = -0.43 [see Fig. [Stb)] .-^ 

To be sure that these changes in Fg and ai are not 
random we carried out systematical studies of the both 
structures at successive decrease of the hole density and 
Drude conductivity. It was recognized that Eq. (|4]) with 
the two fitting parameters, and Fq , describes well the 
experimental data down to CTq — 3.5 ±0.3. All the results 
for ai and Fq are summarized in Fig. 3) The results of 
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FIG. 4: (a) The dependences of Fq''. (b) and (c) The 
(70 dependence of the prefactor ai and the gas parameter r^, 
respectively. Open and sohd circles are the experimental data 
obtained in the present paper for structures 3856 and 3857, 
respectively. The triangles are data from Ref. [2^. The crosses 
are the results of recalculation of the data obtained in Ref. [2ll. 
The curve in panel (a) is theoretical dependence, Eq. The 
curve in panel (b) is the interpolating formula, Eq. ((5]). The 
curves in panel (c) are provided as a guide for the eye. 



Ref. '22' for , obtained for <jq > 30 Go are presented 
in Fig. [DJa) also. One can see that all data match well. 
The scatter of the data from Ref. is broader than that 
obtained here due to the large ballistic contribution that 
complicates the determination of . Note the a vs (Tq 
data can be interpolated by the empirical formula 



a,; = 1 - 



4 Go 

CTO 



(5) 



Let us firstly discuss the behavior of the Fermi liquid 
parameter Fq . Its value as a function of the gas pa- 
rameter, r^, is plotted in Fig. Sfa). It is seen that Fq 
appreciably decreases with that it becomes less than 
—0.454 at Ts ~ 2. It is the value where the interac- 
tion correction in zcroth magnetic field changes sign [see 
Eqs. ([3|) and ([4])]. However, the change of the sign of 
the interaction correction does not result in the metallic- 
like behavior of the total conductivity. This is because 
the insulating-like WL quantum correction dominates in 
our samples. Nevertheless, this fact manifests itself in 
our experiment. Since the triplet channel is suppressed 
with the B-increase, the magnetic field inverts the sign of 
(5(tJ^ again. So the magnetoresistance should be positive 
at low magnetic field and negative at high field. This fact 
graphically shows itself as the maximum in pxx vs B de- 
pendence, which is evident for uo ~ 3.5 Go at B ~ 2.8 T 
[see Fig. [Tfa)]. 

In Fig. SJa) we have plotted the theoretical Ts depen- 
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FIG. 5: (a) The ao dependence of the interaction parameter 
Fq . The symbols are the experimental data. The curves are 
theoretical dependences calculated as described in text. 
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It is seen that experimental points strongly deviate down- 
wards from the theoretical curve with increasing . The 
possible reason of the deviation is the renormalization 
of the Fermi liquid constant F^ with the decreasing 
Drude conductivity which strongly changes with [see 
Fig. mjc)]. This is directly evident from Fig. [51 where 
both the experimental and theoretical [Eq. ^] F^ vs do 
dependences are presented. When calculating the theo- 
retical curves we have used the vs cto dependences from 
Fig. m^c) . It is seen that the lower the Drude conductivity 
the stronger the deviation. 

Theoretically, the effect of renormalization of F^ with 
the changing conductivity was studied in the framework 
of RG theoryj^iii^ii^iii which took the interaction into 
account in the first order in l/cr exactly. According to 
this theory the temperature dependences of both a and 
Fq are the solutions of the system of differential equation 



da 

1^72 



1 + 1 



1 



■72 



72 



ln(l -f 72) 



ill. 
a 



12) 



(7) 
(8) 



where ^ = - In(Tr), 72 = -F^ / (1 + F^), and a is mea- 
sured in units of Go. The term 1 -I- 1 in braces is re- 
sponsible for the weak localization and the interaction in 
singlet channel which in the case of Coulomb interaction 
give equal contributions. 



5 



The above system of differential equations have been 
solved numerically with the following initial conditions. 
We suppose that the high-temperature conductivity is 
equal to the Drude conductivity: (t{^ = 0) = uo. The 
second condition is -f2{£, = 0) = —Fq/{1 + Fq ) where 
Fq is determined by Eq. ^ Note this system describes 
the conductivity as a function of the parameter Tt (i.e., 
as a function of temperature). Experimentally, we are 
able to find the interaction contribution within only the 
relatively narrow temperature range, T = 1.4 — 4.5 K. 
Therefore it is more appropriate to compare the ctq de- 
pendence of Fq rather than the temperature one. The 
solutions obtained for several ctq values, as the Fq vs T 
dependence within the actual temperature range are pre- 
sented in inset in Fig. [5l It is seen that this dependence 
is relatively weak. In order to compare these results with 
the experimental data for Fq we have averaged the calcu- 
lated value of Fq over this temperature interval. Namely 
this averaged value is presented in Fig. [5] as a function of 

CTq. 

It is seen that both the experimental and calculated 
values (labelled as "1st order") of Fq decrease with de- 
creasing ctq- However, the theory gives much faster de- 
crease. Most probably such a discrepancy indicates that 
one should take into account the next terms in 1 /a expan- 
sion in RG equations. To the best of our knowledge this 
has been done only for two particular cases inappropri- 
ate to our situation. The first case relates to multi-valley 
(n^, ^ 1) systems with 72 <C 1.^^ The second one is single 
valley {riv = 1) systems but with the large 72 value^23 

Thus, realizing the crudity of the above estimations we, 
nevertheless, believe that the decrease of the experimen- 
tal value of Fq with the decreasing Drude conductivity 
results from the renormalization of the h-h interaction. 

Strictly speaking, the RG equations ([7]) and ([5]) were 
derived in the absence of magnetic field, whereas the Zee- 
man splitting suppresses the triplet contributions to the 
right-hand side of Eqs. O and ^^t^^i^l Recently, it 
was shown that the effect of Zeeman splitting on conduc- 
tivity can be used for extracting the dependence of Fq on 
temperaturci^i^ii^ In our case the temperature range in 
which the Zeeman splitting is strong, g^sB > T, is small 
fraction of the total interval which we use for averaging. 
Therefore, we do not expect significant difference in our 
results for Fq due to taking into account the Zeeman 



splitting and consider the comparison of our data with 
solutions of Eqs. ([7]) and ^ is almost correct. 

Next we discuss the behavior of the prefactor a^. 
Fig. ^b) shows that decreases sharply when ctq low- 
ers. The behavior of the interaction correction with de- 
creasing (To was studied experimentally for the n-type 2D 
structures in Ref. '2T. The recalculated data from this pa- 
per presented in Fig. Iljb) by crosses demonstrate analo- 
gous decrease also. The possible reason of such vs co 
dependence is the interplay between the interference and 
the interaction which has not been taken into account in 
the RG theory^iii As shown in Refs. [H and [13 two ad- 
ditional terms in the expression for the conductivity arise 
if this interplay is allowed for [see Eq. (40) in Ref. [s^] . 
One term depends on the magnetic field and leads to 
appearance of the prefactor in WL magnetoresistance. 
The second one does not depend on the magnetic field, 
and therefore it was away in Ref. 'st'. It is quite possible 
that namely this term leads to decrease of a^, prefac- 
tor in the interaction correction, with decreasing Drude 
conductivity. Another contribution to the prefactor ai is 
due to the second-loop interaction effect. This correction 
is known for the singlet channel in the unitary ensemble 
(strong magnetic field). '^^ To the best of our knowledge 
the impact of the interplay between the interaction and 
the interference upon the interaction correction to the 
conductivity as well as the second-loop contribution in 
the triplet channel for = 1 is yet to be studied. 

In summary, the behavior of the interaction contribu- 
tion to the conductivity with decreasing Drude conduc- 
tivity is determined both by the renormalization of the 
interaction constant Fq and by the decrease of the pref- 
actor a,; in Eq. (0]), and the latter is more pronounced. 
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